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The Hamiltonian formalism of the Landau-Lifschitz equation for a spin chain with full anisotropy 
is formulated completely, which constructs a stable base for further investigations. 



I. INTRODUCTION 

The Landau-Lifschitz (L-L) equation for a spin chain with full anisotropy is one of the most important completely 
integrable nonlinear evolution equations^. The compatibility pair were first given by E.K.Sklyanin^ in terms of elliptic 
' functions. The Hamiltonian formalism was also examined, but the reduction procedure and results have some open 
— . points so that the final results are not conclusive^. On the other hand, there were some works trying to solve the 
' equation based upon these compatibility pair, but exact solutions were not found until now^. Because these problems 
, are more basic, it is worth to putting a lot of work. In this paper the Hamiltonian formalism is formulated to provide 
' a base for further investigations. 
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II. LANDAU-LIFSCHITZ EQUATION 

The L-L equation for anisotropic spin chain is in the form of 

S f = SxS wx + SxJS, \S\ = 1 (I) 



in which J = diag(Ji, J2, 0/3), J\ < J2 < J3, describing nonlinear spin waves propagating in a direction orthogonal to 
I ■ the anisotropic axis, and the suffices x and t denote the corresponding partial derivatives. Though spin is a quantum 
^ ' quantity satisfying quantum bracket relation, but if one introduces Lie-Poisson bracket for classical spin, all the 
CNj . procedure of disposal is the same as the usual classical onei. The Lie-Poisson bracket is given by 

O ' {S a (x), S b {y)} = -e abc S c (x)5(x - y) (2) 

on : 

where e abc is a fully anti-symmetric tensor, a,b,c= 1,2,3. By using it we write eq.Q in the form of the canonical 
' equation 



d t S a (x)={S a (x),H} (3) 

in which the Hamiltonian is 

/CO 1 
dxH(x), H{x) = ~J2 \(d x S a (x)) 2 - J a S a (x) 2 (4) 
-°° a 

III. JOST SOLUTIONS 

The compatibility pair of the equation are a couple of 2 x 2 matrices given by E.K.SklyaninS, 

L = —iy]w a (\)Saa a (5) 

a 

and 

M = -j ^ U! a (\)SbS cx Craeabc + i Ulb(\)LL>c(\)S a <Ta\eabc\ (6) 

a,6,c a,b,c 



in which 



m 1 m dn(A, k) cn(A, k) 

^l( A )=P— 7T V' W 2(A =p tt r, UJ 3 {\) = p— r (7 

sn(A, K) sn(A, k) sn(A, k) 
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and sn(A, k), etc., are elliptic functions with modulus k: 



J2-J1 



P = t;VJs-Ji (8) 



J3-J1 2 

Since the coefficients ui a (\) are double-periodic functions of the parameter A: 

u) a {\ + AmK + iAnK')=uj a {\) (9) 

where n,m are integers, K and K' are quarter-periods. It is sufficient to consider A in the fundamental period 
parallelogram ft: 

ft : Re A I < 2K, |ImA| < 2K' (10) 

The first compatibility equation is now written as 

d x F{x, A) = L(x, X)F(x, A) (11) 

In the limit of \x\ — > ±00, S — > (0,0, 1), and the asymptotic solution corresponding is E(x, A) = e ~ lu ' 3X ' 73 _ -p^g j og j. 
solutions of eq.(ll) are then defined by 

(jc, A) = ($(jc, A), ^(jc, A)) -► e"'" 3 " 3 , asi^oo 
$(x,A) = (<£(s,A),0(a;,A)) -► e"^ 3 ^ 3 , as a; ->■ -00 (12) 

The monodromy matrix T(A) is given by 

A) = *(*, A)T(A), T(A) = ( "^J j (13) 

V 6(A) a(A) / 



Furthermore, from the periodic properties of sn(A) and sn(A) = sn(A), etc., the compatibility pair (5) and (6) have 
the following reduction transformation properties: 

L{\ + 2K) = a 3 L(X)a 3 , M(X + 2K) = a 3 M(X)a 3 (14) 



L(X + i2K') = a 3 L(X)a 3 , M(A + i2K') = a 3 M(X)a 3 (15) 



IV. LIE-POISSON BRACKET 



From the first compatibility equation, it is found that 

ST(X) 



and 



Defining 5 



SS a (z) 

ST-\X) 
SS b (z) 



= *- 1 (z,X)(iuj a a a )$(z,X) (16) 



= -$- 1 (z,X)(iu b a b )*(z,X) (17) 



{T(X)®T-\X>)} ikijl = {TW&T-Wm} (18) 
the Lie-Poisson bracket of the monodromy matrix is now simply given in the form 

(r W ?r-'(Y) } . e (19) 
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in which the symbol ® in the right hand side is the usual direct product. After substituting eas. l|16|) and (|17|l . the 
explicit expression of eq.(19) is 

where 

R = ^{iwi^o'i <8> i&2 — iu)2w[icr2 ® ci} + Si{ioj2Uj' 3 i(X2 <8> o- 3 — io^a^cs ® &0"2} 

+ S , 2{— W3wi(73 ® CTi + W1W3CT1 ® (T3} (21) 

The eq.(j2UJ has a simple result if the integrand in the right hand side is a full derivative of some function with 
respect to x. In order to do so, we take into account of 

d x (y- l (x, \)a a ${x, A') ®' $- x (x, A')o- a $(x, A)) (22) 

where ctq = /, and a a is Pouli's matrix for a = 1, 2, 3. Here another type of direct product is introduced 

A iro Sy = (A »' B)a Jm (23) 

Using eqs.© and lfTT|) . it is obvious that 

A)$ _1 (a;, A')W a $(x, A)*(x, A') (24) 

in which 

W = i{u a -u' a )S a {a a ®' I -I®' a a ) (25) 

W„ = iS b (uj b a b (7 a + uj' b a a a b ) <g>' tr a + iSWa <8>' (w£(76<7 + uj b a a a b ) (26) 

In fact, ea. H2()(l can be expressed as a linear combination of terms in ea.l|22|l. since R in ea. l)21|l is expressed as a linear 
combination of W a in eas. H25Jl and l|2(i|) : 

i? = /0W0 + /3W3 + /1W1 + / 2 W 2 (27) 

Writing i? and W a in bigger matrices, e.g. 4 x 4— matrices, and comparing the corresponding matrix elements, a 
group of equations for f a are given in the Appendix A. 
Eq.|2n|) is then 

{T(A)f>T- 1 (A')} = ^/aA a (28) 

a 

where 

A a = y\x, A)cr Q *(x, A') ig>' $ _1 (x, A')cr Q $(x, A) ^ (29) 

x— — L 

V. EXPLICIT EXPRESSION OF LIE-POISSON BRACKET 

On the other hand, we have 
3 

/«A tt - /o(A + A3 + Ai + A 2 ) + (/ 3 - / )A 3 + (/1 - / )Aa + (/ 2 - / )A 2 (30) 

Since {T(A)®T _1 (A')} is definite and only the differences fa—fo-, fi~ fo, f2~fo can be determined from eas. (|A4|) ~ (|A6ll . 
we should see that Ao + A3 + Ai + A 2 = 0, which means that the value of fo is of no importance and may be assumed 
to be 0. Thus we have 

{T(A) ® T -i (V)} = /3(A(b) _ A(6Q)) + /l( _ A (&) + A(W)) + / 2 (-A(6) - A(61)) (31) 



in which 



/3 — TTF~. ~T\~i h ~ ~~^7 7\~' h ~ 



2(uj 3 - L0' 3 ) 



2(wi-wi) 



2(W2 — Cl> 2 ) 



and the explicit expressions of A(6), A(60) and A(61) are given in the Appendix B. 
Because of properties shown in (|14() . A can be restricted in the region 

< ReA < 2K, < ImA < 2K' 

In this restriction, A(61) has no contribution, and reduces to 



{T(A)®T- 1 (A')} = / 3 (A(6)-A(60)) 
/ 



Ll>ilu' 2 + U>2W[ 



where ljj, u> ^ m ean u>j (A) , ujj(X'). 
From ca. 1341) . there are 







CJ3— tJg+lO 





— t^g+iO 

i27r<5(w3 - w 3 )|a| 2 






OJ3— W' — ZU 



u 3 — u' — zO 



u>3 
a>3 



{a(A), fe(A')} 



LUiuj'2 + UJ2LO1 



1 



-ab' 



2 — lj 3 + iO 



2(w-a/) 



-ab' 



itt5(uj 3 — u)' 3 )ab' 



MX), b(\')} = - 



; -ab 

2 (^3 — ^3 — 10 

2(w-a/) °" 



177(5(^3 — u>' 3 )ab' 



and then 

{|a(A)| 2 , b(X')} = ~i2TrS{uj 3 ~ u z )uJiw 2 \a\H' 
Furthermore, in the restriction fi_|_, there are 



6{u 3 - u/ 3 ) 



8(X — A') 



as A = A', and 



ujiui'2 + 2 dn(A) dto 3 (X) 



= P 



sn 2 (A) 



= -P- 



dX 



As a result, eqs.(j2SJ), l|3l)|) and (|3Tjl become 



{o(A), b(X')} = ^ + ^ a(A)b(A') + ^(A - A>(A)b(A') 



{5(A), b(A')} = - tt,1 jf a+ T' 5(A)&(A / ) + z^<5(A - A')5(A)b(A') 



and 



{|a(A)| 2 , b(X')} = i2tt P 5(X - X')\a(X)\ 2 b(X') 
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VI. ACTION-ANGLE VARIABLES IN CONTINUOUS SPECTRUM 

As known in inverse scattering transform, a(A) and a(A) are independent of t, and the phase of 6(A) and 6(A) is a 
function of t, which is determined by the asymptotic form of M in eq.@). We have 

b(t,X) = 6(0,A)e- l4wiW2 * (43) 

and then the angle variable is defined as 

Q(A) = arg 6(A) = -In 6(A) (44) 
i 

The action variable P(X) is a function of a(A) and a(A), and usually assumed to be 

P(A)=F(|a(A)| 2 ) (45) 
where F is an unknown function. As these two variables are canonical variables, it should be 

{P(\),Q(\>)} = -5(\-\>) (46) 

By eq.(42), we find 

{F(|a(A)| 2 ), Q(A')} = F'(|a(A)| 2 )27r p( 5(A - A')|a(A)| 2 (47) 
where F' is derivative of F with respect of its argument. Comparing it with eq.(46), we obtain 

F'(KA)| 2 )2 7 rp|a(A)| 2 = -l (48) 



and thus 



Eq.(44) yields 



Hence the Hamiltonian is 



F(A)=F(|a(A)| 2 ) = -^-ln|a(A)| 2 (49) 



Q(A,t) = Q(A,0)-4a; 1 cj 2 (50) 



2K r, n 2K 







pdn(A) 1 
7T J Q ~" sn 2 (A) 



ff= / dA4wi^ 2 P(A) = -- / dA^-^ln|o(A)| 2 (51) 



Therefore, the Hamiltonian has two kinds of expressions: one is an integral with respect to x in eq.Q, and the 
other is an integral with respect to the spectral parameter in ea. i|51|) . Now it is necessary to derive a conservative 
quantity which has two integral forms compatible with the Hamiltonian. 

VII. CONSERVATIVE QUANTITIES 

In the inverse scattering transform, the conservative quantities are derived from the asymptotic form of the first 
compatibility equation. In the limit |A| — > or \k\ — » oo (k = pA _1 ), by using the asymptotic expansion of elliptic 
functions, we have 

L^-i(k + k- 1 r a )S a a a + --- (52) 

and 

fa.ra.ra} = V^{(1 + n 2 ), (1 - 2 K 2 ), (-2 + k 2 )} (53) 
6 

From the first compatibility equation in this limit, writing 

Mx,k)=e lkx+g (54) 
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and expanding 

g x = Vo + (i2fc)~V + (i2k)- 2 r] 2 + ... (55) 

we obtain 

(—iSi + S 2 ) x / 1 c v 
Vo = S 3x + (I-S3) (56) 

— JDi + D 2 

and 

- 2* = ^ + 2 % 2 - (-yi + ^), (??o) + ^ (5?) 

— IDl + D 2 

etcA In general, t/q ^ 0. This situation appears also in the case of isotropic spin chain. In that case, an additional 
phase in the transmission coefficient a(k) was introduced^ to cancel the non- vanishing 770. However, Takhtajan and 
Zakharov pointed out that this is unreasonable 7 . Any way, 771 in ea. (|57|l does not give an expression compatible with 
the Hamiltonian in eq.Q. 

It was formerly shown that the gauge equivalence between the isotropic spin chain and the nonlinear Schrodinger 
equation, which means, by choosing a suitable gauge, the gauge-transformed compatibility pair of isotropic spin chain 
has the same form of that of NLS equation. Therefore, the conservative quantities for the isotropic spin chain can 
be derived from those for the NLS equation by revised gauge transformation and all results expected are naturally 
found 7 -. 

After Takhtajan and ZakharoA* 7 -, it was tried to find some equations that are gauge-equivalent to the L-L equations 
for a spin chain with axial symmetry. However, such equations seem not existent. From a careful analysis of the 
gauge equivalence between the isotropic spin chain and the NLS equation, only the leading terms of the first one of 
compatibility pair are essentials, while other terms corresponding and the second one of compatibility pair are of no 
importance. That is, a gauge is chosen such that it turns the spin in the first order of spectral parameter of the first 
one of compatibility pair into the 3— axis in the spin spaced. 

The explicit expression of the gauge B was given(see Ref.[6]). After the gauge transformation, ea. (|52[) turns to 

L'= -i(k + 2r a S*)a 3 + U + k~ 1 V + ■■■ (58) 

where U = B X B~ X = f ^ q J an d is an 2 x 2 matrix with vanishing diagonal elements. 
Replacing L by L', the similar procedure for eq.JSHJl and (|57Jl yields 

% = 0, v'i = -|«| 2 + 2(r 1 5 1 2 + r 2> S 2 2 +r353 2 ) (59) 

Noticing ©, there are 

{ri,r 2 ,r 3 } = i{J 3 + J2- 2Ji, J 3 - 2J 2 + J u -2J 3 + J x + J 2 } - ~{Ji, ^3} (60) 

since common constant is immaterial 8 . As shown in the case of isotropic spin chain 7 

4M 2 = S ax S ax (61) 

We finally obtain 

^ = i {S ax S ax - {JiSl + J 2 Sl + J 3 Sl)) (62) 
which is just the density of the Hamiltonian in eq.(4). 

VIII. EXPRESSION OF o(A) 

Ea. l59|) indicates a(k) — > as k — ■> 00, so that in this domain 

, ... If „,ln|a(/c')| 2 

lna(fc) = — / dk' / ; 63) 
in J k — k 
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We have 

lna(fc) = / + /i(i2fc)- 1 + --- (64) 

7 =0, h = J dk'\n\a(k')\ 2 (65) 

In terms of A(fc = 2pX~ 1 ), eq.® is re-written in the domain of A ~ 

, If , dn|a(A')| 2 A 

so that 

h = ~ J d\' -^\n\a{\')\ 2 (67) 

It is equal to ea. (|51|l in this domain, since dn(A) ~ 1 and sn(A) aA . 

Since the general formula in this case is invariant with double-periodic transformation, ea. Q66(l may be rewritten as 

lnaMA)) = i/ ^ 3(y) H^3(A0)| 2 (68) 
m J n+ w 3 (A') - w 3 (AJ 

where the integral domain is real axis in f2 + given in eci. (|33[) . We obtain 

lno(w 3 (A)) =I -iI 1 — + ••■ (69) 
2p 



/i = - / ^ 3 (A')ln|a( W3 (A'))| 2 = - / dA'^ln|a^ 3 (A'))| 2 (70) 



where Jo = and 

If d An\a(us(X'))\ 

It approaches to ea. l|51|l in the domain of |A| w 0. 

IX. DISCRETE SPECTRUM 
Ea . (|59|1 should include the discrete part 



where fcj in the complex fc— plane. It can be transformed into A— plane, that is 

A - A,- A 



In general, we write 

^3(A)) = n^^ (73) 

In the limit of |A| w 0, it coincides with ea. (|72|l since U3 (A) pA" 1 in the limit of |A| 0. 
Extending analytically into £!+, the left hand side of ea. (|35fl is 

{lna(w 3 (A)), 6(A fc )} = ■ • • + ^{(ln(w3(A) - u, 3 (A,)) - ln(w 3 (A) - ^(A~))), fe(A fc )} 

3 

= ■ - £ „ ( A)-, 3fe ) {^- 6 < A »» + £ , 3( A)-W ^ b(Xk)) (74) 
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where X k , Xj S fi+ and the right hand side is 

wi(A)w 2 (A fc ) + w 2 (A fe )a;i(A) 1 



2 w 3 (A)-w 3 (A fc ) 

Ea. (|75|l has a pole at w 3 (A) = w 3 (Afc) so that 



b(X k ) (75) 



r c\ \ j,m m wi(A i )w 2 (A fe ) + u 2 (X k )ui(Xj) daj 3 

{wsiXj), b{X k )} = ^ — i — '— ^ — '- — K -^b{X k )5 3k = p— b(X k )5 jk (76) 

Z CLA Xj 

as seen in ea. l|5§|) . Then eq.lO gives 

{X j ,b(X k )} = b(X k )6 jk (77) 

X. ACTION-ANGLE VARIABLES IN DISCRETE SPECTRUM 

We introduce the action-angle variables in discrete spectrum 

Pj = FiuaiXj)), Qj = lnA(w a (Ai)) (78) 



Then we find 



namely, 



dF du3 



(80) 



and then 



du>3(Xj) dXj 

P{w(\j)) = A, (81) 

From eq. 1)73(1 the discrete part of Hamiltonian is 

H d = iA P ^2 (-«s(An) + w 3 (A„)) (82) 



Here A rl lies in f2 + , and the factor 4 stands for four zeros of a(A) in for a single soliton case. We thus obtain 
which is just canonical equation for b(X k ). 



{H d , b{X k )} = iAp—^{Xj, b(X k )} = - ? 4 WlW2 6(A fe ) (83) 

ClXn 



XL CONCLUDING REMARKS 

The Hamiltonian theory to the L-L equation for a spin chain with full anisotropy was examined by E.K.Sklyanin 2 . 
As mentioned, the conservative quantities derived by him are not compatible with the usual conservative quantities, 
such as the Hamiltonian, which affirmatively concludes that his procedure is unreasonable. Moreover, some other 
results may not be beyond doubt. For example, his eq.(2.15) written in the present notation is 

{a(X),b(X')} = -uj 3 (X - X' + iO)a{X)b(X') (84) 

which is questionable, as comparing with eas. (|35|l or H4U|) . 

Mikhailov and Rodin tried to solve the Landau-Lifschitz equation for a spin chain with full anisotropy based upon 
the compatibility pair given by Sklyanin. But explicit solutions were not given. In the isotropic spin case, the Jost 
solutions do not approach to the free Jost solutions as spectral parameter |fc| — > oo. To construct the equations of 
the inverse scattering transform by Cauchy contour integral, Takhtajan introduced a redundant factor to ensure 
the integral having vanishing contribution of the integral along the big circle in complex fc— plane when the radius 
reaches infinity. In the case of spin chain with full anisotropy, the behaviors of the Jost solutions do not approach to 
the free Jost solutions when |A| ~ 0. Rodin and Mikhailov are unable to overcome this difficulty 3 ^. But in order to 
solve the equation, it is necessary to propose a way to do so. 
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APPENDIX A 



To calculate f a in (26), the terms involving S3, Si, S2 in eq.(26) are 



(u)lUl' 2 ai ® (72 — UJ2Ul[<T2 <8 CTl) — 



(h - /o)i(w3 - u' 3 )(I ®' cr 3 - 0-3 I) 

+ (/i - /2)i(^3 + u' 3 )(iai <g>' (J 2 + ia% ®' ai) 



(Al) 



(w 2 CJ 3 (72 ® <X 3 - W3^ 2 0-3 (8 CT 2 ) = (/i - /o)«(^l - ®' <7l - Ol ®' J) 

+ (/2 - /3)i(^i + Wi)(«cr 2 <g>' cr 3 + «0- 3 «' CT 2 ) 



(A2) 



(a^u^cs <8> cri — LJiiv 3 ai ® cr 3 ) 



(/a - /o)i(w2 - c4)( J <8>' cr 2 - (72 ®' I) 

+ {h - /l)*(^2 + W 2 )(zcr 3 ®' (71 + io\ ®' (73) 



(A3) 



among which one equation turns to one another by simple circular permutation of 1, 2, 3. Noticing the direct products 
in two sides are different, and writing them in 4 x 4 matrix form, we obtain 



u>iu; 2 - uj 2 J x = 2(/i - / 2 )(w 3 + ^3), wicj 2 + ci/ 2 cjJ = 2(/ 3 - / )(^3 - 0J3) 



(A4) 



w 2 w 3 - U3uj' 2 = 2(/ 2 - / 3 )(wi + ix4), w 2 w 3 + w 3 w 2 = 2 (/i ~ /o)(wi - o>i) 

W 3 C^ - = 2(/ 3 - /l)(w 2 + W 2 )> w 3^i +^1^3 = 2(/ 2 - / )(W2 - W 2 ) 



(A5) 
(A6) 



APPENDIX B 



Denoting 



A a {L) = m-\L, \)a a V(L, A'), A a (-L) = tf-^-L, A)<r a *(-L, A') 



(Bl) 



We can see 



C a (L) = A> Q $(L, A), C a (-L) = A')(7 a $(-L, A) 



A = A (L) ®' C (L) - A (-L) »' C (-L) = A(6) + A(60) 
A 3 = A 3 (L) (g)' C 3 (L) - A 3 (-L) ®' C a (-L) = A(6) - A(60) 
Ai = Ax{L) ®' Ci(L) - Ai(-L) ®' C7i(-i) - -A(6) + A(61) 
A 2 = A 2 (L) ®' C 2 (L) - A 2 (-L) ®' C 2 (-£) = -A(6) - A(61) 



(B2) 



(B3) 
(B4) 
(B5) 
(B6) 
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where 



Mb) 



( -aV -a'b \ 

-b'a la 1 

-ba' b'a 

V a'b ab' / 



A(60) 



C 



11,11 



ba'e i2( ' UJ:s ~ UJ '^ L b'ae~ i2 ^' 3 ~ u ' 3 " ,L 



a / 6e -i2K-a, 3 )i o C 12 ,21 

V _5' ae i2 ( w 3- w 3) L — 5a' e -*2(w3-"3)£ 



\ 

-afe'e^ 2 ^ 3 "" 3 ^ 
_a'^ e i2 (" 3 -"3)i 

022,22 / 



(B7) 



(B8) 



and 



A(61) = 



(7x1,11 = l'b e -^3^)L _ I h l e i2{u 3 - U ' 3 )L 

Ci2 )2 i = a'ae^ 2 ^ 3 ^ 3 ^ _ a5 / e -i2( W3 - W ^)i 
C21.12 = aW 2 ^-" 3 ^ - aaV 2 ^ 3 -"^ 



22,22 



a'6e i2 ( W3+a, 3)i a 6' e -i2(a) 3 +w 3 )i c U)22 ^ 

l> ae -i2(w 3 +u' 3 )L C*2l,21 -6 a ' e -i2(w 3 +w 3 )i 

C22.ll _o6' e «(u<3+^)i _ a /5 e -i2(a; 3 + W3 )i o y 



(B9) 



(BIO) 



aa e 



/ — i2(w3+w 3 )L 



012,12 = -b'be l2i -^ +u> '^ L + bb'e i2{ " 3+ "* )L 
021,21 = ~b'be- a{ ^ + ^ )L + bb' e - l2{uJ3+ < )L 
022,n = a'ae- j2( " 3+W3)L - aa' e l2{ ^+^ )L 
Here we show the procedure for obtaining eci. (|B3|l as an example. Firstly, there is 

A (L)®'C {L) = 



e -i(w 3 -w 3 )L 



(Bll) 



(B12) 



~/ ae iK-^ 3 )L + lr be -i{^-w 3 )L _gl e i(w' 3 -u 3 )L + y &e -i( U3 -u 3 )L 

+ a'be^ i< - UJ ' 3 ~ UJ;>s,L b'be i ^' 3 ^ U3 ' ,L + a'ae _ '^ 3_W3 ' L 

According to the definition of eg)' in ea. i|23|) . the terms involving vanishing exponent e° and the terms involving 
non-vanishing exponent e ± l2 (^3-^3)£ are collected separately, 



/ a'a -a'b \ 

-b'a b'b 

b'b fc'a 

V a'b a'a/ 



/^ &e -i2(o) 3 -o) 3 )i b'ae- 12 ^-^^ \ 

a / be -i2(u>' 3 -u 3 )L q a / ge -i2(w 3 - W3 )i o 

a'ae 42 ^ 3 "" 3 ^ -a'te ,2 '^- U3 ' L 

V _b' a ^2(u! 3 ^ 3 )L q & /£ e i2(a) 3 - W3 )L y 



(B13) 



Secondly, 



A (-i)®'C (-L) = 



aS / e -i(w 3 -0) 3 )£ + ^/ e i(w 3 -a) 3 )i a ^/ e -i(a) 3 -a) 3 )i _ ^ a / e i(a) 3 -w 3 )L 
fea / e -i( W3 - W3 )£ _ hh i^ 3 -^)L b tf e -i(w 3 -u; 3 )L + ^ e i(w 3 -u' 3 )L 
e -i(ui' 3 -ui 3 )L 







i(w 3 -a; 3 )L 



(B14) 
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is equal to 



( aa! ab' \ / bb'e i2 ^'^ L -l a 'e l2 ^-^ L \ 

66' -la' aa'e- i2 ^ 3 - u '^ L ab'e- l2 ^-^ L 

ba' bV _ a5 / el 2(^ 3 -"3)i go / e i2(u<3-wJ)£ o 

\ -ab' aa' J \ 5a'e~ l2( " 3 ~"3 )L 66' ' e -*2(w3-wJ)i y 



(B15) 



Finally, combining (|B13|) and (|B 1 5|) . we have shown ea. (IB3l 



